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PREFACE 


In  late  1976,  a  study  to  produce  a  wave  climate  for  U.  S.  coastal 
waters  was  initiated  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES).  This  study  was  authorized  "by  the  Office,  Chief  of  Engi¬ 
neers,  U.  S.  Army,  as  a  part  of  the  Field  Data  Collection  Program  which 
is  managed  by  the  U.  S.  Army  Coastal  Engineering  Research  Center.  The 
U.  S.  Army  Engineer  Division,  South  Atlantic*  and  the  U.  S.  Army  Engi¬ 
neer  Division,  New  England,  also  authorized  funds  during  the  initial 
year  of  this  study  (FY  1978)  to  expedite  execution  of  the  Atlantic 
coast  portion  of  this  program. 

This  report  is  the  first  in  a  series  of  four  reports  and  will 
examine  available  data  and  techniques  for  constructing  pressure  fields. 
The  study  was  conducted  in  the  Hydraulics  Laboratory  under  the  direction 
of  Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W. 
Whalin,  Chief  of  the  Wave  Dynamics  Division,  Mr.  C.  E.  Chatham,  Jr. , 
Chief  of  the  Wave  Processes  Branch,  and  Dr.  D.  T.  Resio,  Project  Man¬ 
ager.  This  report  was  prepared  by  Messrs.  W.  D.  Corson,  D.  T.  Resio, 
and  Dr.  C.  L.  Vincent. 

Commanders  and  Directors  of  WES  during  the  conduct  of  the  study 
and  the  preparation  and  publication  of  this  report  were  COL  John  L. 
Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 


SUMMARY 


In  late  1976  a  study  to  produce  a  wave  climate  for  U.  S.  coastal 
waters  was  initiated  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station.  This  climatological  information  is  to  he  produced  by  numerical 
simulation  of  wave  growth,  propagation,  and  decay  under  historical  wind 
fields.  It  is  imperative,  if  such  an  approach  is  to  be  used  for  appli¬ 
cations  of  significant  economic  consequences,  that  the  entire  set  of 
input  data,  all  numerical  techniques,  and  all  general  assumptions  be 
thoroughly  investigated  and  documented  to  determine  the  types  and  mag¬ 
nitudes  of  errors  intrinsic  to  their  use. 

There  are  four  basic  steps  in  the  calculation  of  waves  from  past 
meteorological  data.  First,  pressure  data  must  be  assimilated  into  a 
pressure  field  that  depicts  all  important  synoptic  weather  features. 
Gradients  of  pressure  in  time  and  space,  along  with  certain  thermal 
characteristics  of  the  planetary  boundary  layer,  are  then  used  to  con¬ 
struct  an  estimate  of  a  quasi-geostrophic  wind  speed  and  direction  at 
some  level  where  it  is  assumed  that  the  frictional  effects  of  the  ocean 
surface  on  the  atmosphere  is  negligible.  Next,  an  analysis  of  the  ver¬ 
tical  variation  of  the  wind  in  the  planetary  boundary  layer  is  used  to 
reduce  this  wind  to  a  common  19- 5-m  level.  Finally,  these  surface  winds 
are  input  into  a  numerical  wave  model  to  simulate  wave  generation,  prop¬ 
agation,  and  decay. 

If  any  one  of  the  above  steps  contributes  significant  bias  (on  a 
geographical  basis,  seasonally  or  overall),  it  can  introduce  errors  into 
the  results  that  are  difficult  or  even  impossible  to  remove.  Similarly, 
if  any  step  contains  a  large  random  error,  certain  statistics  (such  as 
duration  curves,  extremes,  and  conditional  probabilities)  can  be  seri¬ 
ously  affected.  Thus,  each  step  must  be  checked  independently  where 
possible.  This  serves  to  substantiate  the  merit  of  the  physics  and  data 
processing  techniques  used  in  each  step  and  hence  tends  to  lend  support 
to  the  worth  of  the  final  product  more  so  than  the  performance  of  only 
wave  comparisons,  regardless  of  how  extensive  these  comparisons  may  be. 
Indeed,  if  each  step  is  shown  to  be  physically  valid,  it  can  be  argued 


that  the  results  should  be  as  accurate  in  sites  where  there  are  no  wave 
data  for  verification  as  they  are  in  areas  where  large  amounts  of  gage 
data  are  available.  Additionally,  if  all  steps  are  modeled  correctly, 
factors  such  as  direction  and  angular  spreading,  which  are  not  generally 
available  for  comparisons,  can  reasonably  be  assumed  to  be  at  least  ap¬ 
proximately  correct. 

There  will  be  four  reports  in  this  series.  The  first  will  examine 
available  data  and  techniques  for  constructing  pressure  fields.  Since 
there  are  few  sets  of  data  available  with  large  samples  of  wind  speeds 
at  the  geostrophic  level,  particularly  since  this  level  fluctuates  in 
accordance  with  the  dynamics  of  winds  in  the  planetary  boundary  layer, 
the  second  report  combines  the  estimation  of  winds  at  the  geostrophic 
level  and  the  reduction  of  these  winds  to  the  surface  level  into  one 
step  which  is  then  compared  with  observations  of  winds  near  the  surface. 
The  third  report  discusses  the  formulation  of  a  wave  hindcast  model  in 
terms  of  nondimensional  growth  characteristics  and  behavior  under  hypo¬ 
thetical  wind  fields.  The  final  report  will  compare  the  hindcast  waves, 
using  the  constructed  oceanic  wind  fields,  with  measured  wave  spectra. 
This  last  report  is  consequently  the  report  that  documents  the  cumula¬ 
tive  effects  of  the  errors  from  all  four  steps  in  the  hindcast  procedure. 
It  is  our  belief  that  numerical  modeling  of  surface  waves  represents  an 
evolution  toward  a  more  reliable  means  of  obtaining  wave  information  for 
climatological  purposes.  Coupled  with  the  concurrent  evolution  of  sta¬ 
tistical  methods,  data  processing  technology,  and  planning  and  design 
capabilities,  this  tool  offers  vastly  improved  ability  to  deal  with 
coastal  problems.  Furthermore,  by  relating  data  to  physical  processes, 
an  underlying  understanding  of  the  wave  phenomena  is  gained.  This  can 
increase  confidence  in  recognizing  the  significance  of  trends,  distri¬ 
butions,  and  correlations  among  various  data  elements,  which  can,  in 
turn,  increase  confidence  in  many  basic  planning,  design,  construction, 
operation,  and  maintenance  decisions. 
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WAVE  INFORMATION  STUDY  FOR  U.  S.  COASTLINES 
SURFACE  PRESSURE  FIELD  RECONSTRUCTION 
FOR  WAVE  HINDCASTING  PURPOSES 


PART  I :  INTRODUCTION 


1.  In  the  derivation  of  wind  fields  from  historical  weather 
records,  it  is  important  to  use  all  information  available  to  the  maximum 
extent  possible,  since  data  density  in  both  time  and  space  can  be  some¬ 
what  sparse.  Simultaneously,  economy  of  effort  is  important  in  a  study 
such  as  this  where  a  20-year  record  of  wind  fields  over  the  Atlantic  and 
Pacific  Oceans  as  well  as  the  Gulf  of  Mexico  is  being  reconstructed. 
Consequently,  some  specialized  sources  of  information  that  are  difficult 
to  incorporate  into  a  machine-oriented  analysis  are  not  considered  as  a 
part  of  this  effort.  For  example,  there  is  a  small  quantity  of  obser¬ 
vations  which  has  not  been  included  either  in  the  initial  construction 
of  the  northern  hemisphere  pressure  fields  or  the  tape  files  containing 
data  from  ships.  This  a  priori  decision  was  made  with  the  understanding 
that  if  it  were  determined  at  a  later  stage  in  the  study  that  signifi¬ 
cant  errors  could  be  avoided  by  supplementing  our  analyses  with  these 
data,  then  these  data  would  be  added  to  the  analyses.  This  report 
describes  the  procedures  used  in  preparing  surface  pressure  field  data 
for  numerical  wave  hindcasting  purposes. 
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PART  II:  AVAILABLE  DATA 


2.  Most  of  the  machine  processible  weather  records  available  for 
the  simulation  of  surface  wind  fields  over  the  sea  come  from  one  of  the 
three  data  files:  (a)  Mil  pressure  fields,  (b)  EDIS  tape  deck  TDF-11, 
and  (c)  EDIS  tape  deck  TDF-lU. 

a.  Mil  pressure  fields.  The  Fleet  Numerical  Weather  Center 
(FNWC)  has  contracted  with  Meteorology  International  Inc. 
(Mil)  to  complete  a  file  of  hemispheric  atmospheric  fields 
on  a  63  by  63  grid  which  overlays  polar  projections  of  the 
Northern  and  Southern  Hemispheres.  Data  obtained  from 
many  sources  are  blended  together  to  form  an  integrated 
file.  Four  observations  per  day  for  a  period  of  more  than 
20  years  are  available  in  this  format.  The  data,  as  pre¬ 
sented,  are  uniformly  distributed  on  the  polar  projection. 

The  observational  data,  on  which  the  gridded  data  are 
based,  are  unequally  distributed  in  time  and  space.  The 
methods  used  for  data  blending  are  described  by  Holl  and 
Mendenhall.* 

b.  EDIS  Tape  Deck  TDF-11.  The  Environmental  Data  Information 
Service  (EDIS )  of  the  National  Oceanic  Atmospheric  Admin¬ 
istration  (NOAA)  collects  shipboard  weather  observations, 
including  air  and  sea  temperature,  wind  speed  and  direc¬ 
tion  and  surface  atmospheric  pressure,  "obtained  from  ship 
logs,  ship  weather  reporting  forms,  published  ship  obser¬ 
vations,  automatic  observing  buoys,  teletype  reports  and 
...  from  several  foreign  meteorological  services."  These 
data  are  sorted  by  Marsden  squares  (ten  degrees  of  latitude 
by  ten  degrees  of  longitude)  and  subsquares  (one  degree  of 
latitude  by  one  degree  of  longitude)  and  by  date  and  recorded 
on  digital  computer  tape.  Data  are  reported  worldwide  but 
are  frequently  sparce  in  many  areas.  There  is  some  redun¬ 
dance  in  the  information  available  from  TDF-11  and  Mil. 
TDF-11,  however,  is  the  only  source  for  the  air  and  sea 
temperature  data  needed  to  define  the  thermal  stability 

of  the  atmosphere. 

c_.  EDIS  tape  deck  TDF-lU.  This  data  file  contains  airway 

surface  observations  including  air  temperature,  wind  speed 
and  direction  and  atmospheric  pressure  adjusted  to  sea 
level.  All  stations  are  on  land  and  only  a  few  are  near 
the  coast,  hence  much  of  the  wind  data  is  not  representa¬ 
tive  of  overwater  conditions. 


*  M.  M.  Holl  and  B.  R.  Mendenhall.  1971.  "The  FIB  Methodology  and  Ap¬ 
plication:  Fields  by  Information  Blending  Sea-Level  Pressure  Version," 
Meteorology  International  Inc.  Report  for  Fleet  Numerical  Weather  Cen¬ 
ter,  Monterey,  Calif.,  under  Contract  No.  N6631j+-70-C-5226. 
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3.  If  an  attempt  were  made  to  estimate  marine  wind  fields  solely 
from  the  wind  reports  contained  in  TDF-11  and/or  TDF-lU ,  the  results 
would  often  be  quite  broadly  interpolated  and  probably  highly  erroneous. 
This  is  due  in  part  to  the  data  sparseness,  in  part  to  the  difficulty 
in  interpolating/extrapolating  a  vector  field,  and  in  part  to  the  over¬ 
reliance  sometimes  on  single  wind-speed  observations  to  fill  in  large 
regions  of  wind  field  over  an  ocean. 

U.  An  alternate  approach  in  obtaining  marine  wind  fields  follows 
from  the  physical  relationships  between  pressure  fields  and  surface  wind 
fields.  On  a  synoptic  time  scale,  the  winds  near  the  earth's  surface 
are  driven  primarily  by  three  factors: 

£i.  The  quasi-geostrophic  wind  located  at  some  level  where 
frictional  effects  of  the  earth's  surface  are  considered 
negligible. 

b.  The  three-dimensional  wind  field  in  the  planetary  boundary 
layer . 

c_.  The  roughness  of  the  surface  over  which  the  wind  is 
flowing. 

Mathematics  that  describe  the  relationships  of  a  parcel  of  air  "to  the 
pressure  gradient,  gravity,  and  the  earth's  rotation"*  can  be  used  with 
an  accurate  representation  of  the  surface  pressure  field  to  calculate  vi¬ 
able  quasi-geostrophic  winds.  The  first  of  the  factors  mentioned  above 
will  be  determined  from  data  developed  by  procedures  described  in  this 
report.**  After  a  preview  of  the  characteristics  of  the  available  data, 
the  surface  pressure  data  from  Mil  were  chosen  as  the  most  applicable  for 
a  starting  point  in  reconstructing  surface  marine  pressure  fields.  A 
detailed  description  of  this  data  set  is  given  by  Holl  and  Mendenhall. t 
5.  In  order  to  determine  if  the  sea-level  pressure  information 
directly  available  from  Mil  was  adequate  for  the  wave  hindcasting 


*  S.  L.  Hess.  1959.  Introduction  to  Theoretical  Meteorology.  Holt, 
Rinehart  and  Winston  Publishing  Co.,  New  York,  p  362. 

**  The  second  and  third  factors  will  be  discussed  in  Report  2  of  this 
series . 

+  See  footnote  on  page  6. 
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purposes  of  this  study,  a  series  of  comparisons  between  Mil  and  National 
Weather  Service  (NWS)  data  were  performed  and  are  presented  in  this  sec¬ 
tion.  The  pressure  fields  for  21  extratropical  storm  systems  (Table  l) 
were  analyzed  to  determine  the  types  and  magnitudes  of  discrepancies  be¬ 
tween  the  Mil  and  NWS  analyses. 

6.  The  storm  systems  were  selected  using  methods  that  should 
not  have  biased  the  comparison  of  pressure  field  representation. 

Storms  1-6  were  the  initial  systems  prepared  for  comparison  and  were 
chosen  from  NWS  synoptic  surface  charts  to  represent  storms  generated 
in  the  Gulf  of  Mexico  or  the  Atlantic  Ocean,  affecting  the  U.  S. 

Atlantic  coast,  and  formed  during  various  seasons.  The  remaining 

15  systems  analyzed  were  selected  as  periods  of  high  winds  and/or 
waves  from  gage  records  for  the  stations  indicated  in  Table  1  and 
shown  in  Figure  1. 

7.  Mil  data  were  available  for  a  63  by  63  grid  covering  a  polar 
projection  of  the  Northern  Hemisphere  (Figure  2)  at  6-  and  2k-hr  in¬ 
tervals.  The  types  of  NWS  surface  charts  used  were  the  3-hr  interval 
North  American  chart  and  the  2U-hr  interval  Northern  Hemisphere  chart. 

8.  The  NWS  Northern  Hemisphere  or  North  American  synoptic  charts 
were  printed  from  microfilm  and  the  Mil  pressure  data  were  contoured  on 
a  drum  plotter.  Since  only  2U-hr  interval  Northern  Hemisphere  synoptic 
charts  were  available  for  the  periods  studied  in  the  ocean  stations 
"India,"  "Juliett,"  and  "Kilo"  areas,  only  2^-hr  intervals  could  be 
analyzed.  Most  other  storms  were  studied  at  6-hr  intervals  (Table  l). 
The  parameters  used  to  examine  the  pressure  data  were  chosen  to  be 
readily  available  from  Mil  sea-level  pressure  data  and  NWS  synoptic 
surface  charts.  After  contour  maps  of  the  Mil  data  and  reproductions 
of  the  NWS  pressure  contours  were  prepared  at  the  same  scale,  the 
storm  geometries  and  pressures  were  studied. 

9.  Storm  geometry  was  compared  by  overlaying  the  Mil  contoured 
pressure  plots  on  the  NWS  charts  (Figures  3-6;  Appendix  B).  As  can 
be  seen  from  Figure  2,  the  Mil  63  by  63  grid  spacing  remains  constant 
over  the  Northern  Hemisphere.  One  grid  space  is  approximately  equal 
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Figure  1.  Location  map  of  data  stations 
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Figure  2.  Mil  grid  on  polar  projection  of  Northern  Hemisphere 

to  381  km  at  60  deg  north  latitude.*  After  visual  comparison  of  storm 
geometry,  it  became  apparent  that  Mil  data  consistently  represented  the 
interior  portions  of  nearly  all,  and  especially  the  more  severe,  storms 
as  less  intense  than  the  NWS  charts  (Figures  3-6;  Appendix  B).  The  grid 
spacing  of  381  km  at  60  deg  north  latitude  is  too  large  to  accurately 
depict  the  interior  portions  of  the  severe  storms  which  are  common  in 
the  higher  latitudes.  Due  to  the  nature  of  the  Mil  grid  and  to  the 


*  Fleet  Numerical  Weather  Center.  197*+. 
Monterey,  Calif. 


"FNWC  Computer  User  Guide," 
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Figure  6.  Mil  and  NWS  pressure  field  overlay  with  cross  section  and 
observed  data  for  Atlantic  coast  Sable  Island  area,  12Z  23AUG73 


typically  less  intense  pressure  gradients  in  lower  latitudes,  the  Mil 
data  representation  of  storms  in  this  region  is  more  similar  to  the  NWS 
representation  than  those  in  more  northern  latitudes  (Appendix  B, 

Plates  B38-B1+2). 

10.  For  Mil  and  NWS  data,  the  storms'  lowest  millibar  contours 
were  recorded  as  the  central  pressures  (CP).  Plots  of  the  NWS  CP  versus 
Mil  CP  show  that  Mil  consistently  represents  the  storm  central  pressures 
higher  than  the  NWS  (Appendix  C).  During  analysis  of  the  central  pres¬ 
sures,  the  observed  sea-level  pressure  data  presented  on  the  NWS  syn¬ 
optic  charts  were  compared  with  Mil  and  NWS  pressure  field  presentations 
(Figures  3-6).  When  possible,  the  history  of  the  more  severe  storms  was 
reviewed  in  the  Mariners  Weather  Log  and  the  Monthly  Weather  Review  to 
discern  if  later  analysis  had  changed  the  pressure  field  presented  by 
the  NWS.  As  seen  in  Figures  3-6,  the  NWS  sea-level  pressure  data  were 
consistently  found  to  be  compatible  with  observed  data  and  hence  are 
accepted  as  valid  representations  of  storm  pressure  fields.  On  the 
other  hand.  Mil  discrepancies  from  the  NWS  pressure  fields  were  found 
to  be  unsubstantiated  and  thus  must  be  considered  as  a  source  of  error 
in  any  geostrophic  wind  speeds  calculated  from  these  pressure  fields. 
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PART  III:  DEVELOPMENT  OF  PRESSURE  FIELDS  FOR  HINDCASTS 


11.  Once  it  was  determined  that  the  Mil  pressure  fields  did  not 
suffice  for  a  synoptic-scale  representation  of  pressure  gradients  (and 
consequently  not  for  an  adequate  wind  field  description  for  hindcast¬ 
ing),  the  question  arose  as  to  how  to  compensate  for  this  problem  in  a 
manner  that  was  both  accurate  and  economical.  After  considerable  effort 
in  reviewing  storms*  in  the  Atlantic  Ocean,  it  was  determined  that  most 
of  the  problems  with  the  Mil  pressure  fields  occurred  in  the  western 
Atlantic  near  the  east  coast  of  the  U.  S.  and  Canada.  The  reason  for 
this  is  that  this  region  is  one  of  rapid  cyclogenesis,  storm  develop¬ 
ment,  and  storm  movement;  therefore,  storm  features  tend  to  be  smaller, 
more  variable,  and  more  changeable  in  time  here  than  toward  the  eastern 
Atlantic.  Since  there  did  not  appear  to  be  a  reliable  means  of  trans¬ 
forming  the  Mil  pressure  field  into  a  NWS  pressure  field  in  this  region, 
it  was  decided  that  the  NWS  pressure  field  would  be  digitized  and  that 
this  digitized  information  would  be  overlaid  on  the  Mil  data,  with  some 
blending  on  the  edges  to  achieve  some  degree  of  continuity  between  the 
two  fields.  Figure  7  outlines  the  steps  taken  to  produce  the  pressure 
field. 

12.  The  base  grid  used  in  pressure  field  development  is  the  area 
bounded  by  Mil  coordinates  I  (28  to  63)  and  J  (l  to  32)  on  the  63  by  63 
Northern  Hemisphere  grid  (approximately  the  southeast  quadrant  shown  in 
Figure  2).  This  grid  area  covers  all  of  the  North  Atlantic  and  portions 
of  the  bordering  continents. 

13.  As  was  shown  in  Part  II,  the  spacing  of  the  Mil  grid  (381  km 
at  60  deg  north  latitude)  is  inadequate  for  the  representation  of  impor¬ 
tant  synoptic-scale  storm  features.  As  a  first  step  in  increasing  the 


*  The  term  storm  is  used  in  this  report  when  referring  to  extratropical 
cyclones  (atmospheric  depressions  or  low  pressure  systems)  typically 
located  along  a  front  that  separates  two  different  air  masses.  No  NWS 
pressure  data  from  tropical  cyclones  (hurricanes  or  typhoons)  are  in¬ 
put  into  the  pressure  field.  These  storms  will  be  treated  in  a  sepa¬ 
rate  analysis  for  wave  heights  and  are  thus  excluded  from  entering 
into  the  pressure  analysis  here. 


REVIEW  WEATHER  CHARTS- 


Mn  PRESSURE 


DIGITIZE 
> COASTAL 
REGION 


RESPECIFY 
CENTRAL  PRESSURE 
FOR  OCEAN  STORMS 


*-  EDIT,  REVIEW 


PRODUCE 

FINE 

GRID 


PRESSURE  GRID 
COLLATION  PROGRAM 


Figure  7-  Pressure  field  development  procedure 

resolution  of  the  pressure  field,  the  32  by  3 6  base  grid  is  subdivided 
by  a  factor  of  U  in  both  directions,  creating  a  125  by  lHl  grid  into 
which  additional  information  is  blended  (Figure  8). 

I1*.  Figure  9  shows  that  portion  of  the  Atlantic  in  which  storms 
appeared  to  have  characteristics  more  appropriately  represented  by  the 
NWS  analyses  (shown  as  Coastal  Grid  Border).  This  is  also  the  region 
in  the  Atlantic  in  which  storms  affect  the  coastal  U.  S.  in  a  local 
sense,  i.e.,  through  sea  and  not  through  swell.  Consequently,  this 
region  is  chosen  as  a  fixed  area  that  will  be  digitized  if  necessary. 

It  is  not  necessary  to  digitize  every  weather  map  and  blend  this  in, 
however,  since  not  every  day  contains  storm  features  of  the  type  that 
are  not  well  represented  by  the  Mil  gridded  pressures. 

15-  As  noted  previously,  the  Mil  pressure  data  become  increas¬ 
ingly  inaccurate  as  smaller  scales  are  approached.  Along  the  U.  S. 
Atlantic  coast,  many  small  storms  produce  waves  that  can  be  significant 


Mil  GRID  INCREMENTS  381  KM 


WIS  GRID  INCREMENT  =95  KM 

Figure  8.  Mil  grid  divided  into  spacing 

in  terms  of  a  climatological  representation  of  coastal  processes.  In 
an  effort  to  ensure  an  adequate  description  of  these  smaller  storms  as 
well  as  to  achieve  a  better  description  of  the  larger  storms,  the  fol¬ 
lowing  criterion  was  chosen  as  an  objective  guide  to  the  selection  of 
weather  maps  to  be  digitized. 

16.  The  fixed  area  shown  in  Figure  9  would  be  digitized  when  a 
storm  consisting  of  at  least  two  closed  isobars  (at  H-mb  intervals)  was 
located  with  a  center  within  this  area.  All  storms  of  even  moderate 
intensity  that  are  overly  smoothed  by  the  Mil  analysis  should  be  re¬ 
solved  in  this  process. 

17-  The  fixed-area  grid  is  aligned  at  1+5  deg  to  the  base  grid 
in  order  to  cover  the  area  desired  and  to  ease  interpolation  to  a  finer- 
than-MII  grid  spacing  (Figures  9  and  10).  NWS  pressure  data  within  the 
coastal  grid  are  digitized  from  North  American  synoptic  surface  charts 
into  a  33  by  21  array.  However,  this  array  does  not  have  data  points 
at  all  1/1*  Mil  grid  locations.  The  array  is  interpolated  to  a  65  by 
1*1  array  that  furnishes  the  data  locations  required  to  transfer  the 
NWS  data  into  the  1/1+  Mil  grid  system  (Figure  10).  Only  NWS  data  are 
used  within  the  fixed-area  grid.  Interpolation  to  Mil  data  occurs  only 
at  the  grid  borders  where  the  pressure  data  are  linearly  interpolated 


LONGITUDE  LATITUDE 


Figure  9.  Coastal  grid  on  polar  projection 


over  one  Mil  grid  space  from  the  fixed-area  grid. 

18.  Digitization  of  the  NWS  pressure  field  provides  a  good  rep¬ 
resentation  of  coastal  storms  in  the  Atlantic  Ocean;  however,  storms 
outside  of  this  area  can  be  significantly  weakened  in  the  Mil  represen¬ 
tation  as  compared  with  the  NWS  representation.  In  order  to  intensify 
storms  in  the  Mil  representation  in  a  simple  manner  to  be  consistent 
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oOZ 


COASTAL  GRID 


Figure  10.  Coastal  grid  with  subdivisions 


with  the  NWS  representation.  Northern  Hemisphere  synoptic  charts  are  re¬ 
viewed  for  major  storms  not  within  the  coastal  grid.  These  storms  are 
selected  (or  not  selected)  according  to  the  storms'  intensity  inferred 
from  central  pressure  recordings  at  various  latitudes  (Table  2).  The 
selection  process  is  designed  to  provide  a  simple,  consistent  method  for 
selecting  storms  that  may  not  have  been  represented  by  Mil  as  severe  as 
they  should  have  been.  This  selection  method  should  allow  for  the  re¬ 
specification  of  the  more  severe  storms  and  some  of  the  less  severe 
storms.  Crossplots  of  the  selected  Mil  CP  versus  the  corresponding  NWS 
CP  for  the  first  year  analyzed  show  that  most  storms  selected  during  the 
review  process  required  intensification  (Figure  ll),  which  tends  to  vali¬ 
date  the  selection  methodology. 

19.  Once  the  storms  to  be  respecified  in  the  Atlantic  (outside  of 
the  coastal  grid)  have  been  selected,  the  central  pressure  from  these 
storms  is  blended  into  the  Mil  pressure  field  representation.  This  is 
accomplished  by  developing  a  Ho  by  Uo  subgrid  (in  the  expanded  1/k  Mil 
scale)  about  the  location  of  the  Mil  minimum  pressure  (CPMII)  and  sym¬ 
metrically  blending  the  NWS  central  pressure  into  the  Mil  grid  using  the 
following  algorithm: 


PB  = 


X1(CPNWS} 


X2Ps 


where 

X1  =  (d2/r)2 

X2  =  (d^r)2 

P„  =  blended  pressure  for  a  given  location 
a 

CPnwS  =  central  pressure  reported  by  the  NWS 

Pg  =  pressure  at  a  given  location  in  the  1/k  Mil  grid  system 

d^  =  distance  from  the  boundary  radius  where  the  deepening 
factor  goes  to  zero  (5  Mil  grid  units)  to  the  point 

d2  =  distance  from  the  location  of  the  lowest  Mil  pressure  (the 
assumed  center  of  the  low)  to  the  point 

r  =  sum  of  d^  +  &2  which  must  equal  the  total  radius  over 
which  blendig  is  performed 


Since  d2  goes  to  zero  at  the  boundary  radius  the  pressures  become  equal 


.  _ ■  . 
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1000 


to  the  Mil  pressures  there,  and  since  d^  goes  to  zero  at  the  assumed 
center  of  the  low,  the  pressures  become  equal  to  the  NWS  central  pressure 
there;  in  between,  a  distance-squared  interpolation  is  affected. 

20.  Figure  12  presents  a  sample  of  xnitial  and  blended  pressure 
fields.  The  blended  pressure  field  produces  a  steeper  pressure  gradient 
than  that  developed  by  the  Mil  pressure  field.  Since  the  blended  pres¬ 
sure  field  is  dependent  on  the  Mil  data  for  placement  and  field  geometry, 
some  misplacement  (with  respect  to  the  NWS  placement)  of  the  blended 
pressure  field  does  occur. 
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Example  of  storm  central  pressure  respecification 


PART  IV:  CONCLUSIONS 


21.  Based  on  analysis  presented  in  PART  II  of  this  report,  the 
sea-level  pressure  field  developed  by  the  NWS  is  accepted  as  a  standard 
with  which  other  data  should  be  compatible.  The  pressure  field  developed 
from  Mil  data  is  not  usually  as  intense  as  the  NWS  pressure  field  in 
areas  with  steep  pressure  gradients.  The  grid  system  in  which  Mil  data 
are  computed  and  stored  does  not  appear  to  be  able  to  depict  the  interior 
portions  of  severe  storms.  In  order  to  develop  a  valid  pressure  field 
for  use  in  the  Atlantic  Coast  Wave  Information  Study,  the  Mil  pressure 
data  are  supplemented  (when  required)  with  NWS  pressure  data  using  two 
types  of  subgrid  input  systems:  a  fixed-area  or  coastal  grid  and  a 
center  respecified  grid.  Since  our  investigations  show  that  the  NWS 
pressure  charts  are  consistently  in  agreement  with  historical  observa¬ 
tions,  the  resulting  surface  pressure  fields  from  this  procedure  should 
constitute  an  excellent  starting  point  for  hindcasting  oceanic  wave 
climates. 


2h 


'd 

o 

s* 

cd 


in 

<D 

p 

in 

£ 


■3 

Eh 


cd 

cj 

•rH 


CM 


c3 

a> 


cd  cd 
<d  a> 


crt 

crt 

crt 

fl 

crt 

crt 

5 

■3 

>3 

cj 

a> 

a; 

a> 

a; 

<L> 

<L> 

<L> 

a> 

> 

u 

u 

u 

H 

u 

P 

p 

P 

crt 

crt 

crt 

H 

o 

■3 

< 

■3 

< 

*3 

■3 

a 

c 

d 

crt 

QJ 

QJ 

<u 

3 

H 

0 

d 

0 

<l) 

in 

u 

• 

• 

• 

• 

• 

• 

o 

o 

O 

H 

■3 

3 

3 

£ 

• 

in 

cn 

in 

cn 

cn 

in 

in 

cn 

£3 

e 

fl 

3 

P 

cn 

P 

p 

M 

M 

M 

M 

w 

M 

3 

3 

o 

£ 

£ 

£ 

P 

M 

cn 

cn 

(U 

<u 

<u 

£ 

a 

cd 

cd 

<u 

crt 

crt 

<D 

<u 

0 

a> 

<1> 

a> 

OJ 

OJ 

OJ 

O 

•H 

•H 

•H 

•H 

cn 

o 

o 

rH 

H 

rH 

rH 

rH 

rH 

rH 

xs 

*d 

^d 

H 

d 

CJ> 

o 

P 

P 

P 

p 

p 

P 

p 

p 

p 

•H 

d 

d 

d 

0 

tui 

crt 

crt 

crt 

crt 

crt 

crt 

p 

p 

p 

H 

M 

M 

*“D 

o 

o 

CO 

CO 

CO 

uj 

OJ 

OJ 

o 

o 

o 

£ 

£ 

£ 

£ 

£ 

3 

•H 

•rH 

o 

o 

o 

P 

P 

r 

* 

r 

r> 

r» 

• 

• 

• 

• 

• 

a 

-p 

P 

P 

c 

p 

a 

P 

P 

p 

P 

P 

P 

n 

* 

r\ 

CO 

CO 

CO 

CO 

CO 

P 

cn 

cn 

cn 

aJ 

cn 

crt 

cn 

cn 

cn 

cn 

cn 

in 

p 

p 

p 

10 

od 

cd 

crt 

crt 

m 

crt 

crt 

crt 

crt 

crt 

in 

in 

in 

• 

• 

• 

• 

• 

al 

o 

o 

o 

P 

o 

P 

o 

o 

o 

0 

o 

0 

d 

crt 

cd 

o 

o 

o 

o 

o 

O 

o 

o 

o 

c 

o 

<: 

o 

o 

O 

o 

u 

o 

o 

o 

o 

O 

o 

o 

o 

r> 

rs 

* 

u 

o 

o 

nd 

o 

tJ 

cj 

o 

cj 

o 

o 

a 

CJ 

o 

o 

o 

o 

o 

•H 

•H 

•H 

c 

•H 

d 

•H 

•rH 

•rH 

•H 

•H 

a 

a 

CJ 

•H 

•rH 

•rH 

•H 

•H 

•H 

P 

P 

P 

cd 

P 

crt 

P 

p 

P 

P 

P 

P 

•H 

•H 

•H 

P 

P 

P 

P 

p 

P 

d 

d 

d 

d 

a 

d 

d 

d 

d 

d 

<W 

Ch 

d 

c 

d 

d 

d 

d 

cd 

cd 

aJ 

<H 

cd 

Cm 

cd 

aJ 

a} 

crt 

crt 

crt 

•H 

•rH 

•H 

crt 

crt 

crt 

crt 

crt 

cd 

H 

H 

rH 

rH 

rH 

H 

rH 

rH 

rH 

rH 

rH 

rH 

O 

CJ 

O 

rH 

rH 

H 

rH 

rH 

rH 

P 

p 

P 

d 

P 

0 

P 

P 

P 

P 

P 

P 

crt 

crt 

crt 

P 

P 

p 

P 

P 

P 

< 

<c 

< 

o 

< 

o 

< 

< 

c 

< 

<e 

c 

(i. 

p-, 

< 

< 

<: 

C 

< 

<d 

a 

. 

. 

• 

• 

• 

. 

• 

• 

• 

• 

• 

• 

• 

• 

p 

p 

p 

p 

p 

. 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

P 

P 

P 

p 

P 

CO 

d 

d 

?H 

u 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

O 

O 

O 

O 

O 

• 

ID 

tD 

ID 

ED 

u 

u 

u 

u 

u 

u 

ID 

£5 

£5 

ID 

a 

cn 

a 

0 

o 

•H 

p 

■P 

cd 

r 

> 

!h 

0) 

f» 

in 

u 

P 

4) 

O 

P 

C 

M 

CM  CM  CM  CM  CM  CM  CM  CM  CM  CM 


O 

o 

O  vo 

o 

o  vo 

CM 

CM 

CM  VO 

o 

O 

o 

CM 

CM 

CM 

CM 

CM  CM 

o 

o 

o  o 

o 

o 

o 

rH 

rH 

«H 

O 

o 

o 

o 

rH 

rH 

H 

H 

H  rH 

rH 

o 

ir\  c*— 

CM 

ON  CM 

~=t 

rH  OO 

t/N  t— 

rH 

1 — f 

CM 

00 

LA 

LA 

LA  t'- 

cn 

on 

H 

H  O 

rH 

rH 

CM 

CM 

o 

CM 

CM 

H 

CM 

rH 

CM 

rH 

CM 

rH 

H  CM 

v  * 

o 

rH 

CM 

OJ  CO  CO 

o 

O 

CM 

»H 

rH 

p'- 

00 

LA 

rH 

CM 

On  rH 

p 

fr*. 

rH 

H 

o  o 

o 

O 

o 

O 

O 

rH 

rH 

rH 

H 

O 

o 

O 

O 

rH 

o 

O  H 

crt 

3 

PO 

on 

On  iH 

o 

00 

p- 

00 

CO 

oo  on 

OO  CM 

00 

00  VO 

vo  vo 

CM 

•H  H 

Q 

3 

VO  VO  vo  C—p- 

t- 

VO 

t— 

t— 

t —  t—  t— 

t — 

t— 

LTV 

LA 

LA  VO  VO  VO 

« 

1 

1 

1  1 

rH 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i  i 

o 

o 

o 

O  CM 

<H 

o 

CM 

o 

CM 

O  CO 

o 

CM 

o 

o 

CM 

CM 

CM 

CM 

OJ  MD 

o 

o 

O  H 

O 

o 

rH 

o 

rH 

O 

H 

O 

H 

o 

o 

rH 

rH 

rH 

rH 

H  O 

o 

s 

OS  C— -st  VO 

00 

rH 

C—  CM 

OO 

i — I 

C— -st 

lA 

o 

o 

0\  On  On  Ox  H  VO 

p 

5h 

CM 

o 

H  O 

c- 

rH 

rH 

CM 

CM 

o 

CM 

CM 

H 

CM 

rH 

rH 

O 

rH 

O 

H  CM 

CO 

o 

1 — 1 

r*--sr 

OJ 

CM  ao  CO 

o 

O 

CM 

rH 

rH 

00 

00 

LA  rH 

CM 

On  H 

rH 

rH 

o  o 

o 

O 

O 

O 

rH 

H 

H 

•H 

O 

o 

o 

O 

rH 

O 

O  H 

00 

00  ON  H 

00 

-=r 

oo 

00 

oo 

00 

00  CM 

00 

oo  -4-  vo  vo 

CM 

rH  H 

vo  \o  vo  r— 

c- 

VO 

t—  c— 

tr— 

c— 

LA 

LA  LA  VO  VO  VO 

CM  00  -St  LA  VO  CO  ON  O  H  CM  CO  -st  lAVD  t-CO  OnO  H 
HHHHHHrHHHHCMOi 


Year /month/ day/hour . 


Table  2 

Selection  Guidelines  for  Storms  Not  Within  Coastal  Grid 


For  North  Latitude 

Central  Pressure 

20°  to  40° 

<990  mb 

1+0°  to  50° 

<980  mb 

50°  to  60° 

<970  mb 

60°  to  65° 

<960  mb 

Note:  No  storms  located  over  land  or 
latitude  are  selected. 

above  65 0  north 

APPENDIX  A:  SYSTEM  HISTORY  FOR  SELECTED  STORMS 


Source:  National  Weather  Service  North  American  Synoptic  Surface  Charts 


Time 

GMT 

Day 

CP 

mb 

Location 

Lat ,  Long . 

Comments 

18-31  October  1963 

1200 

18 

1012 

25°N,73°W 

Attached  to  a  warm  front 

0000 

19 

1008 

27°N,72°W 

1800 

20 

992 

32°N,75°W 

0000 

21 

9&k 

3l+°N,75°W 

l+-mb  contour  spacing  about  0.25°  of 
latitude 

0600 

21 

981+ 

3l+°N,75°W 

Storm  assigned  as  Hurricane  "Ginny" 

1200 

28 

1000 

3l+°N,75°W 

l+-mb  contour  spacing  about  2°  of  latitude; 
still  named  "Ginny" 

1200 

29 

981+ 

1+1°N,67°W 

"Ginny"  attached  to  a  cold  front;  50+  knot 
wind  speeds  reported 

1800 

29 

972 

i+i+°n,6o°w 

and 

36°N,70°W 

Storm  developed  double  low;  no  longer 
designated  "Ginny" 

0000 

31 

988 

Southern  low  moved  inland;  contour  spacing 
about  1°  of  latitude 

6-10 

November  1963 

0000 

6 

1008 

33°N,80°W 

1800 

6 

992 

35°N,77°W 

Area  from  981+  contour  to  1012  contour 

covered  U.  S.  Atlantic  coast 

1200 

7 

981+ 

37°N,77°W 

0600 

10 

996 

l+l°N,62°W 

12- 

■15  July  1969 

0000 

12 

1012 

36°N,75°W 

1200 

13 

1000 

1»20N,68°W 

0000 

15 

1001+ 

1+3°N,60°W 

5-7  April  1971 

1800 

5 

1012 

28°n,81+°w 

0600 

6 

1008 

32°N,79°W 

A2 


Time 

GMT 

Day 

CP 

mb 

Location 
Lat .Long. 

Comments 

5-7  April 

1971 

(Continued) 

-• 

1200 

6 

1000 

34°N,78°W 

4-mb  contour  spacing  from  1000  mb  to 

1012  mb  was  about  1°  of  latitude 

0600 

7 

992 

39°N,72°W 

13-14  January  1973 

0000 

13 

1012 

Storm  elongate  (NE-SW)  off  southeast 

U.  S.  coast 

0600 

13 

1004 

30°N,75°W 

1800 

13 

996 

31°N,71°W 

25-  to  30-knot  wind  speeds  reported  along 
North  Carolina  coast 

0600 

Ik 

1000 

34°N,64°W 

30-knot  wind  speeds  reported  on  Bermuda 

9-12  February  1973 

1200 

9 

1008 

26oN,90°W 

1800 

9 

1004 

27°N,87°W 

0000 

10 

1004 

Storm  moved  across  Florida;  40-knot  wind 
speeds  reported  off  North  Carolina  coast 

1200 

10 

992 

32°N,77°W 

50-  to  75-knot  wind  speeds  were  reported 
along  North  Carolina  coast 

1800 

10 

988 

32°N,77°W 

0600 

12 

992 

39°N,67°W 
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APPENDIX  B:  COMPARISONS  OF  Mil  PRESSURE  FIELDS  AND  NWS  PRESSURE  FIELDS 

1.  The  contour  interval  for  Plates  Bl-BlU  is  4  mb. 

2.  The  contour  interval  for  Plates  B15-B37  is  5  mb. 

3.  The  contour  interval  for  Plates  B38-B42  is  U  mb. 
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APPENDIX  C:  CROSSPLOTS  OF  Mil  CENTRAL  PRESSURES 
AND  NWS  CENTRAL  PRESSURES 
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PLATE  Cl 


NWS  CP  VERSUS  Mil  CP 
FOR  STORMS  1  -  12,  &  21 


In  accordance  with  letter  fron  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Corson,  William  D 

Wave  information  study  for  U.  S.  coastlines;  Report  1:  Sur¬ 
face  pressure  field  reconstruction  for  wave  hindcasting 
purposes  /  by  William  D.  Corson,  Donald  T.  Resio,  Charles  L. 
Vincent.  Vicksburg,  Miss.  :  U.  S.  Waterways  Experiment  Sta¬ 
tion  ;  Springfield,  Va.  :  available  from  National  Technical 
Information  Service,  1980. 

24.  [7]  p.,  [23]  leaves  of  plates  :  ill.  ;  27  cm.  (Tech¬ 
nical  report  -  U.  S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  ;  HL-80-11,  Report  1) 

Prepared  for  Office,  Chief  of  Engineers,  U.  S.  Army, 

Washington,  D.  C. 

1.  Atmospheric  pressure.  2.  Climatological  data.  3.  Mathe¬ 
matical  models.  4.  Meteorological  data.  5.  Numerical  simu¬ 
lation.  6.  Water  wave  hindcasting.  7.  Water  waves. 

I.  Resio,  Donald  T. ,  joint  author.  II.  Vincent,  Charles 
Linwood,  joint  author.  III.  United  States.  Army.  Corps  of 
Engineers.  IV.  Series:  United  States.  Waterways  Experiment 
Station,  Vicksburg,  Miss.  Technical  report  ;  HL-80-11,  Report  1. 
TA7.W34  no. HL-80-11  Report  1  _ 


